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Abstract 

Context. Photoelectric heating is a dominant heating mechanism for many phases of the interstellar medium. We study this mecha- 
nism throughout the Large Magellanic Cloud. 

Aims. We aim to quantify the importance of the [C II] cooling line and the photoelectric heating process of various environments in 
the LMC and to investigate which parameters control the extent of photoelectric heating. 

Methods. We use the BICE [CII] map and the Spitzer/SAGE infrared maps. We examine the spatial variations in the efficiency of 
photoelectric heating: photoelectric heating rate over power absorbed by grains, i.e. the observed [C II] line strength over the inte- 
grated infrared emission. We correlate the photoelectric heating efficiency and the emission from various dust constituents and study 
the variations as a function of Ha emission, dust temperatures, and the total infrared luminosity. The observed variations are inter- 
preted in a theoretical framework. From this we estimate radiation field, gas temperature, and electron density. 
Results. We find systematic variations in photoelectric efficiency. The highest efficiencies are found in the diffuse medium, while 
the lowest coincide with bright star-forming regions (~1.4 times lower). The [CII] line emission constitutes 1.32% of the far in- 
frared luminosity across the whole of the LMC. We find correlations between the [C II] emission and ratios of the mid infrared and 
far infrared bands, which comprise various dust constituents. The correlations are interpreted in light of the spatial variations of the 
dust abundance and by the local environmental conditions that affect the dust emission properties. As a function of the total in- 
frared surface brightness, S'tir, the [C II] surface brightness can be described as: S[cii] = 1-25 S^m [10 erg s^^cm~^sr~^], for 
Stir -c^ 3.2 ■ 10~*erg s~^cm~'^sr~^. We provide a simple model of the photoelectric efficiency as a function of the total infrared. 
We find a power-law relation between radiation field and electron density, consistent with other studies. The [C II] emission is well- 
correlation with the 8 /im emission, suggesting that the polycyclic aromatic hydrocarbons play a dominant role in the photoelectric 
heating process. 

Key words. Galaxies: Magellanic Clouds - ISM: dust, extinction - Galaxies: infrared - Galaxies: lines and bands 



1. Introduction 

The structure and evolution of the interstellar medium (ISIVI) 
is largely dependent upon the therrnal processes taking place 
([Goldsmith et al."1969^,'de Jond^mt iMcKee & Osrika^flOTTt 
iDra ine 1978; Ferriere et al. 1988). This, in turn, shapes the evo- 
lution of galaxies as a whole, as the constituents of the ISM 
are responsible for the characteristics of incipient stellar genera- 
tions. Therefore, an understanding of the agents which dominate 
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the heating and cooling of interstellar gas is of fundamental im- 
portance. 

A dominant heating source of the ISM of galaxies is the pho- 
toelectric (PE) emission of interstellar dust grains. Absorption 
of a far-ultraviolet (FUV) photon by a dust grain may result 
in the ejection of an energetic electron which heats interstellar 
gas via collisions. Photoelectric emission as a he a ting m ech- 
anism of the ISM was first proposed by Spitzer (1948) and 
later revisited by Watson ( 1972) and de Jong ( 1977). Since then, 
it has been found that the process dominates the heating of a 
range of interstellar media: neutral atomic gas clouds, the photo- 
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diss ociation regions (PDRs), and the warm inter-cloud me dium 
(e.g. lMaciel & Pottaschlll982l:IWeingartner& Drainell2001h 

The PE heating proces s has received much theoret- 
Watsonl ^ 



TgTl jdeJond 11977 



iDraind 
1994t 



ical attention 

1978 ^, 'Tielens & Hollenbach' '1985b' 'Bakes & Tiekni 
Dw ek & Smith 1996; Weingartner & Draine 2001). Due to grain 
charging, PE heating efficiency is highly dependent on the phys- 
ical conditions which determine the ionisation and recombina- 
tion rates. Specifically, it depends on the FUV radiation field, 
gas temperature and electron density. In turn, the extent of grain 
charging and therefore the efficiency of PE heating is also highly 
dependent on the grain species involved. 

The C"'' fine structure transition CP 3/2 

the dominant coolant of the diffuse ionised 

atomic gas as well as i n PDRs (Ipalgarno & McCra 



-'P1/2) is 
and diffuse 



19721: iTielens & Hollenbaciil Il985blal: IStacey et a .1 I199U 
Madden et all 119931: iPetuchowski & Bennett! 119931: iHeilesI 



1994 '). Several reasons account for its dominance: carbon is the 
fourth most abundant element and it has an ionisation potential 
of 1 1.3 eV, less than that of H. This C+ transition is also easy to 
excite as it has a relatively low excitation temperature (~92K). 
Thus, it is able to cool war m neutral gas (T = 30 — lO^'lC) 
whereas other species can not ( Tielens & Holl enbach 1985b a; 
IWolfire et al. 1990). The efficiency of the C+ as a coolant is de- 
pendent upon environment. When temperatures or densities are 
high, other li nes, primarily [O I] 63jj,m, par tic ipate in the cool- 
ing p rocess dTielens & HoUenbachl Il985bllat [Hollenbach et alj 
Il99lh . The critical density of the transition is relatively 
low (~ 3 • lO'^ cm^"^). At densities above the critical density or 
temperature above 92 K the cooling by the [C II] line saturates 
and its importance as coolant diminishes. 

An observational study of PE heating and gas cooling re- 
quires [CIlQ and infrared (IR) observations, covering wave- 
lengths tracing the variety of dust components in the ISM, and 
ideally of sufficient spatial resolution to separate environments, 
because of the dependence of the process on environment and 
composition. Because of its proximity, the Large Magellanic 
Cloud (LMC) is an obvious candidate for a study of these pro- 
cesses. We undertake such a study using the Spitzer legacy pro- 
gram: Surveying the Agents of a G alaxy's Evolution (SAGE 
iMeixner et al] 120061: [Bernard et al] 1008), a nd the Balloon- 
borne Infrared Carbon Explorer mission (BICE iMochizuki et alj 
I1994 . SAGE fully mapped the LMC at high spatial resolu- 
tion from 3 to 160/im while the latter offers a [CII] map of 
the entire LMC. These datasets offer advantages over previ- 
ous work because of their enhanced spatial resolution, wave- 
length coverage, and sensitivity. Prior studies using NASA's 
Kuiper Airborne Observatory (KAO) and the Long Wavelength 
Spectrometre (LWS) aboard the Infrared Space Observatory 
(ISO) e xamined [C II] emission integrated acro s s whole galax- 
ies (e .g. IStacevetalJll991l: iMadden et alj|1993l: llVIalhotra et all 
2001 ), or in specific region s within galaxies (e.g. Stutzki et alJ 
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19961: IMadden et al.l ll997^. There have not been many [C II] 



studies which probe the range of different phases of the ISM. 
Moreover, even though PE efficiency is believed to be a strong 
function of grain size, there have not been many prior observa- 
tional studies of the effect of PE heating due to distinct grain 
populations. Most studies have considered the total PE heating 



Tiu'oughiout tliis paper, we utilise tlie following notation: we refer 
to the fine structure transition with and to the emission line that it 
produces with [C II]. 



across all grain populations because they did not have access to 
bands which trace the emission from distinct grain populations. 

The composition of the ISM of the LMC makes it 
an interesting la boratory because it has a low met alhcity; 
Z ^ 0.3 - 0.5Zq IWesterlundl (1 19971) and Z ^ O.25Z0 iDufoui! 
(119841) . The presence of metals, in the form of dust, is integral 
to the PE effect. How does the low metallicity lower dust abun- 
dance affect PE heating? 

Moreover, it is known from observational studies that the 
LMC and other low metallicity galaxies have a dearth of poly- 
cyclic aromatic hydrocar bons (PAHs) compared to Galactic 
value s (ISakon et al.l 120061: [V ermeii et al. 2002b; Madd en et al 
20061 IWu et all 120061: iGalliano et al.. ,2007; .Engelbracht et al 



20081) . This condition raises another question: since it is found 
theoretically that PAHs should play the greatest role in the PE 
heating process, how is the PE heating affected in a galaxy with 
a prominent dearth of PAHs? 

In contrast to the metallicity argument given above, it is 
found from observations that the C+ coolant plays an even 
more important role at low metallicity than it does at high 
metallicity , contributi ng up to ~10 times more to the far- 
infrared (FIR) emiss ion dPoghtsch et al.lll995l:llsrael et al.ll 19961: 
Ma dden et al.lll997 ^. The larger relative strength has been ex- 
plained by the clumpy nature of the ISM in low metallicity 
galaxies. Observations of [C II]/CO show that for low metallic- 
ity galaxies, this ratio can be 10-30 times higher than for normal 
galaxies or active galaxies (e.g. Madden 2000). The far UV radi- 
ation penetrates deeper into the molecular cloud at low Z, for the 
same Ay, leaving a smaller CO core and larger C+ emitting en- 
velopes. The result is a preponderance of CO cores clumps with 
larger exposed surface area. As the C"*" resides primarily in the 
envelopes of the clouds the increase in surface ar ea results in a 
higher ratio of [C II]/CO emission (e.g. lRoUig et aLl(l2'007l) ). 

The aim of this paper is to explore the qualitative behaviour 
and observationally quantify the extent of PE heating and [C II] 
cooling in relation to environment. The format of this paper is as 
follows: Sec.|2]reviews the observational details of the data used 
in this study. Section |3] presents data treatment process and the 
final images. 

Section |4| quantifies the importance of the [C II] coolant 
globally across the LMC. In Sec.|5]we examine PE heating as 
a function of environment based on the Ha surface brightness 
and Sec. |6]explores the variations in [CII] emission within the 
LMC using the Ha criterion and dust temperature. To explore 
the correlation between [C II] emission and the emission from 
the various grain components in the LMC, Sec. [Tj compares the 
spatial distribution of the [C II] emission and the emission at 
various Spitzer bands. The key-concepts of PE heating are re- 
viewed in Sec. 17.21 which are applied in Sec. |8]to describe the 
PE heating as a function of radiation field within the LMC. Using 
the observed relation we calculate electron densities, and find a 
correlation between radiation field and electron density (Sec.|9]l. 
The paper concludes with Sec. [10] an analysis of the dependence 
of PE heating on grain population; we study the qualitative be- 
haviour of the extent of PE heating on grain population, and we 
also quantify the contributions to PE heating from the various 
grain populations. 



2. Observations 

Relevant observational details of all the maps used in this study 
(the [C II], Spitzer and Ha maps) are summarised in Table [1] 
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6^ 5^30"^ 5^ 4^30^" 

Right Ascension {2000) 

Figure 1. (a) and (b) images convolved and re-gridded to the resolution and pixel scheme of the [C II] map. To enhance 

the contrast of the low level emission, each image has been truncated at 0.25 times the max of that particular image. The [C II] 
map is contoured over the images (thin grey and white contours) at levels of 1(T[cii]. The highest contours are white so that they 
can be distinguished from the underlying images. The thick grey lines enclose the regions which are defined as star forming given 
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6^ 5^30'^ 5*^ 4^30 

Right Ascension (2000) 



Figure 2. 70/im (a) and 160/im (b) images convolved and re-gridded to the resolution and pixel scheme of the [CII] map. To 
enhance the contrast of the low level emission, the images have been truncated at 0.25 times the max of that particular value. The 
[C II] map is contoured over the images (thin grey and white contours) at levels of Icrjc n] • The highest contours are white so that 
they can be distinguished from the underlying images. The thick grey lines enclose the regions which are defined as star forming 
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Line/Band 


[C II] 


8 nm 


24 fim 


70 nm 


160 jim 


Ha 


Mission 


BICE 


SAGE 


SAGE 


SAGE 


SAGE 


O T T A A 

SHASSA 


la limit [MJysi--^] 




0.2 


0.2 


1.0 


2.0 




la limit [IQ-^erg s-^cm^^sr^M " 


0.47 


3.0 


0.7 


1.0 


0.8 


1.21 • 10-2 


Ict limit [MJysr-^] " 




4 • 10-* 


4 ■ 10-* 


4 ■ 10-3 


2 • 10-2 




la limit [IQ-'^erg s^^cm-^sr-^" 


0.47 


6- 10-3 


1 ■ 10-3 


5 ■ 10-3 


1 • 10^2 


4.45 • lO-'^ 


Survey Area 


6° ■ 10° 


7.1°- 7.1° 


7.8°- 7.8° 


7.8°- 7.8° 


7.8°- 7.8° 


13° • 13° 


Beamsize ["] 


894 


2 


6 


18 


40 


240 


Linear size [pc] 


~225 


~0.5 


~1.5 


~4.5 


~10 


~60 


Reference 


[1] 


[2] 


[2] 


[2] 


[2] 


[3] 



(1200 Ih 



2.1. The [CW] line 



2.3. Ha Emission 



The entire LJVIC was mapped in the 158/im [CII] Hue by 
iJMochizuki et al.l (Il994j) during the BICE mission. The veloc- 
ity integrated [C II] line surface brightness map was continuum 
subtracted by using a linear baseline and foreground subtracted 
using COBE data to estimate the MiUcy Way contribution. The 
measurements are calibrated against observations of M17 carried 
out bv Matsuhara et al. (1989) . We compare several region s of 
the 30 Doradus co mplex mapped bv lPoglitsch et alj (Il995h and 
llsrael et al.l (Il996h aboard the KAO withe the BICE map and find 
agreement better than ^ 20%. This is within the 30% calibra- 
tion un certainty of the BICE map as quoted bv iMochizuki et al] 
(11994. The one a lev el for the BICE map as determined by 
llVIochizuki et alj (Il994l) is 0.47 • lO^^erg s^^cm^^sr^^ 

The BICE map offers a 6° x 10° field of view of the LMC 
(see also Figs.[T]and|2l). The beam has a FWHM of 14.9', corre- 
sponding to a linear size of ^ 225pc at the distance of the LMC 
(50 kpc Feast 1999; Keller & Wood 2006). The positional un- 
certainty is ~ 6' ( Mochizuki et aL,1994.) . We have regrided the 
original map into pixels of 3' 45pc) in length. 

A significant fraction of the pixels are noise dominated: 
46.1%, 21.7% and 10.8% of all pixels in the final map are above 
are above 1, 2 and 3 Sa^c nj, respectively. The contours in Figs.[T] 
and |2] are spaced at multiples of (T[q uj . To sample the [C II] 
emission from the diffuse ISM, in our analysis throughout the 
rest of this paper, we consider pixels above the 2(7[c ii] level. 



2.2. Infrared Dust Emission 

To probe dust properties and abundances, we use the SAGE 
infra red and mid-infrared survey of the LMC (iMeixner et al.l 
I2OO6 ). The survey offers a 7° x 7° field of view of the LMC 
at eff ective wavelengths of 3.6, 4.5, 5.8 8.0 (IRAC iFazio et all 
12004 . 24.0, 70.0 and 160.0Aim (MIPS iRieke et alJl2004 . We 
use the 8, 24, 70 and 160/im bands as the most important dust 
emission components can be traced primarily by these bands. 



Zodiacal and Milky Way contamination light at MIPS wave- 
lengths has been removed by doing off-source subtraction. The 
IRAC data have not been background subtracted as the back- 
ground contribution is negligible compared to emission from the 
LMC. Averaging several regions located off the LMC in the 8/im 
image yields a background contribution of about O.G4MJysr^^. 
Such a level is indeed small compared to the mean value of about 
0.46MJysr^^, for the 8/im image. 



As a general way to distinguish physical environments within 
the LMC , we use an Ha surface brightness criterion. The Ha 
map used is part of the Southern H-alpha Sky Survey Atlas 
(Gaustad et al. 2001). The LMC is imaged in a single 13° x 13° 
field and does not require mosaicing. The images are continuum 
subtracted and smoothed to a resolution of 4'. The sensitivity 
level of the Ha maps is about 0.5 Rayleigh ( = 1.21 • 10^^ erg 
s^^ cm'^ sr"^ at Ha). 



3. Data Treatment and Presentation 

All IRAC, MIPS and Ha data were convolved to the lowest res- 
olution data, the [C II] map. The shape of [C II] beam is not 
precisely determined, and we therefore convolved the data with 
a Gaussian k ernel with a FW HM of 14.9', the FWHM of the 
BICE beam ( Mochizuki et al.al994.) . After convolution, the data 
were interpolated to match the [C II] pixel scheme. 

To determine a values for the maps after the data treatment, 
we carried out a facsimile of the data treatment process with sim- 
ulated Gaussian noise images. Standard deviations of the noise 
images after the full data treatment were measured, and are in- 
cluded in Table[T] We use those a values in our analysis. 

Maps of the four Spitzer bands (convolved and re-gridded) 
with overlays of the [C II] contours are presented in Figs. [1] 
and |2] It is evident from Figs. [T^ and ^ that the 8/im and 
160/im emission extends significantly into the bulk of the galaxy. 
Figs.[T}5 and[T^ show that the 24/im and 70/im emission is much 
more concentrated toward H II regions such as 30 Dor and Nil. 
More attention to the spatial distributions of the grain emission 
will be given in SecQ 

Note the slight offsets of the [CII] emission peaks from 
the IRAC and MIPS emission peaks (especially near Nil), 
which can be larger th an the 6' poin t ing ac curacy. Similar de- 
viations were noted by Kim & ReachI (l2002h in their analysis of 
atomic gas in conjunction with the BICE map. We have carefully 
checked for problems in the coordinate encoding used in the var- 
ious maps. The observed displacements do not appear to be due 
to such problems and thus we conclude that there are real offsets 
between the peaks of emission of [C II] and IR emission, which 
can be as much as 100 pc. Indeed, the region near Nl 1 is known 
to have most of its molecular material in the direction of the dis- 
placement as compared to the illuminating sources. Perhaps the 
observed displacement is due to offsets in the peak emission at 
the molecular interface, as such displacements a re known to also 
occur in the Galaxy (e.g. lCesarskv et an(ll996t) '). 
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4. [C II] Emission Giobaliy Across thie LMC 



5. Distinction of Phiysicai Environments 



For comparison with other galaxies, we calculate the ratio of 
[C II] to the total infrar ed (TIR) integrated over the LMC. We 
adopt the expression of iDale & Helou (2002) to calculate the 
TIR from the Spitzer the 24, 70 and 160/^m filters, approxi- 
mating the integrated 3 to 1000/im infrared surface brightness. 
Stir- 



Stir = 1.559zyS,.^24 + 0.76 86z/S 



1^,70 



160- 



(1) 



The coefficients were derived from measured i nfrared SEP 
shapes of a sample of galaxies observed by Spitzer. [Praine & Lil 
(|2007|) also provide an equat ion for the TIR lumin osity using a 
modified prescription of the iDale & Heloul (|2002|) SED model 
and incorporating the Efim IRAC band as well as the MIPS 
bands. The difference between the two methods is less than 15%. 
Since both prescriptions were made using SED properties from 
entire galaxies, it is important to check that they yield good esti- 
mates on smaller scales as well. We therefore compare the results 
from Eq. ([T]i with a straight integration over all the IRAC and 
MIPS bands. We find no systematic differences in the value of 
TIR as derived from the two methods for distinct environments, 
i.e. SF and the diffuse medium. Moreover, the values agree to 
within ~10%, which can most likely be accounted for by the 
differences in the precise wavelength ranges that both calcula- 
tions consider. 

The Spitzer definition of TIR emerged from the definition of 
LpiR which was motivated by the IRAS bands (Helou 1986) and 
covers the FIR wavelength range from 42 to 122/im. . The differ- 
en ce between TIR and FIR has been observationally quantified 
bv lHunter et al.l(l200ll) who find that for irregular galaxies (such 
as the LMC), Ltir/Lfir ~ 2, (with a dispersion of only a few 
percent). 

Integrated across the entire galaxy, we find that the to- 
tal [CII] luminosity in the LMC is (5.9 ± 1.8) ■ lO^Lp as- 
suming the distance to the LMC to be 50kpc (iFeaslI 11995 
iKell er & Wood" 20060 consistent with iMochizuki et al ' 



andJCim & Reach (2002)who estimate Ljcuj in the LMC to be 
5.7 • IO'^Lq and 6.5 • IO^Lq respectively. We find that the value 
of L[cn]/LTiR is 0.46 ±0.14%. Assuming a factor of 2 between 
TIR and FIR, L[cn]/LFiR = 0.9%. 

The relative contribution of [C II] to the integrated FIR 
has often been used to evaluate the globa l star formation 
activity in galaxies (e.g. IStacev et al.l Il99lh . This value for 
the LMC is high compared to normal and gas rich galax- 
ies which normally have values of L[cii]/Lfir less than 
1%. Values of L [ r;TT]/LFiR = ■}% to 0.5 % are typical 
(IStacev et alJll99ll; iMalhotr a et al."l997', '2001); for the Milky 
Way L[cii]/Lfir ~0.3% (Wright etal. 1991). Low metallic- 
ity galaxies can typically have L[r^TT]/LFiR as high as 1% to 



roicaiiy 

3% (iPoditsc h et al.lll995l:llsrael et alii996tlMadden et al.ll997t 
iMad den 

200^" This higher ratio is a consequence of the low 
metallicity: due to the reduced dust abundance, the overall mean 
free path of UV photons can be larger, resulting in a decrease in 
the FIR intensity arriving at the surfaces of the molecular clouds. 
To add to this effect, the lower dust attenuation results in the C+- 
emitting regions being larger as the photo-dissociating photons 
traverse a larger volume of the molecular cloud. 

^ The contribution of the [C II] emission to the infrared emission is 
usually quoted in ratios of luminosities and we therefore adhere to this 
convention. We note, though, that we actually calculate values for this 
ratio with surface brightness. This is, however, equivalent to a ratio of 
luminosities assuming that both the [C II] and TIR are radiated isotrop- 
ically. 



We aim to study the PE heating and the [C II] cooling line as 
a function of environment. While the spatial resolution of 225 
pc results in some mixing of phases, we can still delineate dis- 
tinct average conditions. We define two environments using the 
Ha surface-brightness. Ha emission is efficient at distinguishing 
between dense H II regions and diffuse media. This is because 
H II regions are Ha bright due to their high free-electron densi- 
ties, which results in a high recombination rate. The Ha line is a 
tracer of the physical parameters which determine the degree of 
ionisation. That is, radiation field, temperature and density; all 
of which play critical roles in PE heating and the C+ transition. 

iKim & ReachI (|2002|) studied the LMC using the BICE 
[CII] map and made a similar distinction to separate the 
phases of the ISM. They studied the C+ cooling rate for re- 
gions with an Ha surface brightness, Sho, above and below 
4.25 • lO ^^erg s~^cm~^sr~- '^. This number was based on the 
work of iKennicutt & Hodgd (Il986l) . We considered using this 
criterion, but this causes several well known H II regions such 
as N41, N144 and N132 to be classified as diffuse. We conclude 
that the criterion should be lowered to better represent the dif- 
ferent phases of the IS M. For a reformulation of th e criterion, 
we examined the data of IKennicutt & Hodgd (Il986l) . They pho- 
tometrically observed HII regions in the LMC and tabulated 
the Ha surface brightness. A histogram of the distribution of 
these surface brightnesses, shows the rapid fall off at the low- 
est values, suggesting that the faint side of the distribution is 
noise-dominate d, while the bright s i de is comprised of reliable 
values. Indeed. IKennicutt & Hodgd (11986') warn that their mea- 
surements of H II regions with the lowest surface brightnesses 
are unreliable. We therefore take the value of the peak of the dis- 
tribution, 1 ■ lO^^erg s~^cm~^sr~^ , as the lowest reliable sur- 
face brightness for an H II region, which we use to distinguish 
between physical environments in the LMC. Every pixel with an 
Ha surface brightness below this value we call "diffuse", and 
every pixel with an Ha surface brightness above it, we call "star 
forming" (SF). 

The thick grey lines in Figs. [T] and |2] enclose the SF pixels, 
while the diffuse pixels reside outside the grey lines. One can see 
from this figure, that the SF regions correspond to the brightest 
HII regions, such as 30 Dor, Nil and the prominent HII re- 
gions along the LMC bar, while not extending too far into the 
diffuse medium. While we can not avoid including diffuse emis- 
sion within these SF regions, most SF pixels are dominated by 
H II regions and PDRs. Likewise the regions we label diffuse 
will undoubtedly contain some denser ionised material, but will 
be dominated by the diffuse conditions. 

To validate the threshold, we estimate the electron density 
(rie) for the sets of SF and diffuse pixels under the Case B ap- 
proximation. The approximation provides Shq which depends 
upon electro n density, size of t he emitting regions and gas tem- 
perature (e.g. lValls-G abaud'1998). For the estimation, we use the 
mean values of Sho for the SF and diffuse pixels, and a tempera- 
ture of lO'^K. The sizes of the emitting regions were determined 
by examining the original Ha images and visually determining 
the physical sizes of typical SF regions and the voids between 
them. For the SF pixels, we assume that the light is dominated 
by the emission from H II regions. The lengths determined were 
about 5' and 15' for the SF and diffuse regions respectively. For 
the SF pixels, we assume that the light is dominated by the emis- 
sion from H II regions. This yields an average iio of ~ 100 and 
~ lcm~'^ for the SF and diffuse regions, respectively. The latter 
density is reasonable for densities of the warm ionised medium 
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( Nordgren et alJ[T992l) . The former is consistent with lPeck et al.l 
(Il997h . who find a mean electron density of ~ 200 cm~^ in 
30 Dor. 



6. Variation of [C II] Emission Withiin thie LMC 

To examine the variations of L[c nj /Ltir, we show a ratio map 
of L[cii]/Ltir overlaid with the boundaries of the SF regions 
in Fig. [3] The pixels with the lowest L[cii]/Ltir are associated 
with the centres of the brightest SF regions, confirming the sce- 
nario in which the PE heating is least efficient at high density. 

The distribution of [C II] and TIR from the different phases 
are summarised in Table |2] Although the SF regions have much 
higher values of [C II] surface brightnesses, approximately half 
of the [C II] emission originates from the diffuse medium. The 
SF regions contribute just less half of the [C II luminosity 
(~48%). 

Fig. |4] shows L[cii]/Ltir as a function of the Ha surface 
brightness. As can be seen, the efficiency is roughly constant 



across the range of Ha surface-brightnesses. There are some de- 
viant pixels with high Ljcuj/Ltir, in particular at the lowest 
values of 5hq, which are most-likely the result of the noise in the 
[C II] map. The bulk of our points clutter around Ljdtt] /Ltir — 
0.005, whic h is similar t o the v alues found by iJVIalhotra et al.l 
(1200 ll) and iHunter et all (1200 lb . who examine spiral and ir- 
regular galaxies (including the LMC). There is evidence that 
on the whole the 30 Dor region has a slightly lower ratio of 
L[cii]/Ltir, than the rest of the LMC by less than 10% (see 
also Fig.O. 

Theoretically, such a decrease in this ratio towards 
the densest regions is exp e cted. The PDR m odels of 
iTielens & Hollenbachl(fT9"85blal) : IWolfire et al.l(ll99a) show that 
[C II] emission levels off at the highest gas temperatures, radia- 
tion fields and gas densities. Other lines take over (part of) the 
cooling process at high gas temperature and high gas density. 
Since dust temperature roughly scales with gas temperature and 
gas density, it is possible that the observed decrease is associ- 
ated with the critical temperature and critical density of the C+ 
transition being reached. We, however, do examine alternative 
explanations in the following sections. 

7. Distribution of [C II] Emission and Grain 
Component Emission in thie LMC 

The 8, 24, 70 and 160/im bands used in this study trace the 
emission from distinct grain populations. We use these measure- 
ments to study [CII] cooling and PE heating with respect to 
grain abundance. We adopt the generally accepted interpretation 
of the dust constituents and at which w avelength they emit (e.g. 
iDesert et al.lll990l; iDraine & Lill200ll) . The 8^m band is domi- 
nated by PAH emission. The 24/Ltm band mainly traces the emis- 
sion of stochastically heated carbonaceous grains with sizes less 
than ~ 0.01/im; termed very small grains (VSGs). The 160/im 
band mainly traces the grains larger than ^ 0.01/im, in radia- 
tive equilibrium; termed big grains (BGs). Finally, the emission 
detected in the 70/im band probably represents a combination 
of the BG continua and the VSG emission, and may also trace 
grains stochastically heated, thus not in thermal equilibrium with 
the interstellar radiation field. Hereafter, we represent the PAH, 
VSG and BG components with 8, 24 and 160 /im band, respec- 
tively. 

The spatial variation of PAH emissi on in galaxies 
has been studied extensive ly (e.g. Roche & Aitk enI 119851: 
"Xitken & Roche 1985; Lea chI 1 1 9871; IVoitll 19921; ICes arskv et^aD 
1996; Verstraete et al.. .200lt ISiebenmorgen et al.. .20041 
Povich et al. 2007). The studies consistently find a lack of PAH 
emission for the most active regions, i.e., H II regions, starburst 
galaxies and AGNs. An interpretation is that the PAHs in these 
regions are destroyed due to the hard, intense radiation field. 
The VSG grain emission peaks in the H II region while the PAH 
emission peaks in the adjacent PDR. 

We present ratio maps of Sg, S24, S70 and Sigo to Stir in 
Fig. |5^-d. The grey lines enclose the SF regions as given by 
our Ha criterion. Consistent with the studies mentioned above. 
Figs. |5^,b show that the relative PAH emission is weak in the 
SF regions, while the VSG emission peaks in the SF regions. 
Fig-EJl shows that the BG emission relative to TIR mostly peaks 
in the diffuse medium. BGs may be present throughout the LMC, 
the high flux of UV photons in the H II regions entails that they 
are hotter in these regions and therefore radiate at wavelengths 
blueward of the 160/im band. Finally, we note that the 70/im 
emission relative to TIR (Fig. |5};) most closely follows that of 
the VSG emission. 
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Table 2. Contribution of [C II] from various regions in the LMC. 



Region 


L[CII]/L[C II], total 


Area/Areatotai 


LtIr/LtIR, total 


L[c ii]/Li6o 


L[cii]/LfTiR] 


L[cii]/LpiR] 


Diffuse regions 


51.8 ±0.1% 


86% 


59. 0± 0.1% 


0.032±0.1% 


0.40±0.01% 


0.80±0.01% 


SF regions 


48.2 ±0.1% 


14% 


41.0 ±0.1% 


0.063±0.1% 


0.55±0.01% 


1.10±0.01% 


30 Dor 


8.1 ±0.1% 


0.8% 


9.0 ±0.1% 


0.072±0.1% 


0.42±0.01% 


0.84±0.01% 


Nil 


2.2 ±0.1% 


0.7% 


2.0 ±0.1% 


0.058±0.1% 


0.57±0.01% 


1.14±0.03% 


total LMC 


100% 


100% 


100% 


0.042±0.1% 


0.46±0.01% 


0.92±0.02% 



" Also referred to as the photoelectri c efficiency (see Sec . 7.2 ). ''FIR is defined using the IRAS bands (see sectior [Z2l l. TIR assumes the factor of 
2 increase from the FIR, as found bv lHunter et alj ( 1200 Ih and discussed in the text (section |2j2]l. Due to the calibration uncertainties of the BICE 
map, all of these values have a systematic uncertainty of ±30%. The uncertainties in the table reflect the random noise which is important for 
comparing the relative. 
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Figures. Ratio maps of the 8 (a), 24 (b), 70 (c) and 160/im (d) bands to the TIR emission. Contours enclose the SF regions. (See 
online version for colour.) 



7. 1. Grain Component Emission Relative to [C\\] 
Emission 



To study how the [C II] emission varies with grain population, 
we examine this parameter as a function of the ratios of grain 
component emission. It should be noted, though, that the ratios 
represent ratios of grain emission and not that of grain abun- 
dances. 

Fig.|6]presents S[cii] as a function of S8/S24, S24/S160, and 
Ss/Sieo. We omit discussion of ratios with the 70 /im band as 



they exhibit the same behaviours as the ratios with the 24 ^m 
band. We note the following trends: 

- In the SF regions where the [C II] emission is highest, the 
PAH emission is low relative to the VSG emission (Fig.|6h), 
which is explained by the fact that the PAHs do not survive 
in H II regions where the VSGs are known to peak. As the 
S8/S24 ratio increases, the pixels trace more and more of the 
diffuse medium where the [C II] emission is the lowest. 

- Panel b shows again that those regions with very prominent 
24 fim emission, i.e. SF regions, show strong [C II] emission. 
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Figure 6. [C II] surface brightness (top) and the S[c ii]/Stir (bottom) as a function of grain component emission ratios as traced by 
the Spitzer bands. The black dots and red crosses denote the different diffuse and SF pixels, respectively. Note, that the high values 
in the bottom three panels are due to noise in the [C II] map (See text for details). See online version for colour. 



- Fig. |6}; shows that there is no correlation between the [C II] 
emission and the Ss/Sigo ratio for the diffuse medium. Some 
fraction of the SF pixels, though, show a rough increase with 
[CII] as a function of Sg/Sigo- The regions with the most 
intense radiation fields will heat the surrounding grains to 
very hot temperatures so that the BG emission will shift to 
bluer wavelengths, outside of the 160/im filter. Toward the 
SF regions, the amount of PAH emission also decreases (see 
above). The decrease in the 160/im band is greater than the 
decrease in the PAH emission toward the most intense SF re- 
gions. This distinction between the brightest and the fainter 
SF regions is not reflected in the relative [C II] strength 
(Fig. IS). 



7.2. Theoretical description of PE 

PE heating occurs as follows: absorption of a far-ultraviolet 
(FUV) photon by an interstellar dust grain liberates an elec- 
tron within the grain. The electron travels through the grain, es- 
capes, and then overcomes any Coulomb attraction if the grain 
is charged. If the absorbed energy exceeds the work function of 
the grain plus its Coulomb potential, the electron escapes with 
excess kinetic energy. That energy goes into heating the ISM via 
collisions with the gas species. Because of this Coulomb poten- 
tial, PE heating efficiency is highly dependent upon the charge 
state of the grain. This, in turn, is dependent upon the physical 
conditions which determine grain ionisation and recombination 
rates. The PE efficiency (e) is defined as: 



where Pabs, represents the power absorbed by the grains and 
FpE is the photoelectric heating rate. Thus, e is the fraction of 
the power absorbed by the grains that goes into heating the ISM. 

Photoelectric efficiencies as a function of environmental 
conditions, grain size distributions and grain c ompositions were 
calculated semi-empirically by de Jon g|(ll977), and a b initio ca\- 
cul ations have been ma de by iWeingartner & Draind (12001 1) and 
Bak es & TielensI (|1994|) . 

Assuming an MRN g rain distribution and only carbonaceous 
grains. B akes & Tielen i (fT994,) derive an analytic expression 
for PE efficiency (valid for gas temperatures much less than 
10 000 K): 



e(Go,T,nc) 



4.87-10-2 

1 + 4 • 10-3^0.73 



(3) 



PE 

abs 



(2) 



where 7 is the ratio of ionisation to recombination rate 
(GoT^/^Hc^^). Thus, high ionisation rates will increase grain 
charge and therefore lower efficiency. 

7.3. A Proxy fore 

L[cii]/Ltir is a proxy for e given several assumptions: i) that 
PE heating and the 0+ transition dominate the heating and cool- 
ing processes respectively, i.e. FpE = A = ^[cii] (where A rep- 
resents the cooling rate) and ii) And that interstellar grains and 
molecules re-radiate all of the energy absorbed in the infrared. 
However, the [CII] line is not always the dominant coolant. 
Therefore, to observationally estimate e, some a uthors include 
the measured luminosities of other FIR lines (e.g.'M eixner et all 
[1992; Young Owl et al. 2002). We do not have maps of the LMC 
in other FIR lines, and thus can not include their contribution in 
the cooling rate. Therefore, it should be kept in mind that our 



10 



Rubin et al: [C II] and dust emission in the LMC 



calculations of e represent a lower limit on the actual values of 
PE efficiency, in particular in the densest regions. 

Here we estimate the contribution from other lines that may 
be missing within the 15' beam. We perform this estimation 
in the 30 Dor region, as this is the region in the LMC where 
the contributions of other FIR lines to the cooling rate should 
be the most significant. We estimate the contribution from the 
[OI] 63/im line, the dominant cooling line in regions where 
[C II] cooling is suppressed. The 3 Dor complex was mea- 
sured in [C II] and [O I] (63nm) bv iPoghtsch et al.1 (Il995b with 
a 55" (FWHM) beam aboard the KAO. They found peak in- 
tensities in [C II] and [O I] of 1 • 10~'^ erK s~^cm~^sr~^ and 
6.1 • lO^^erg s^^cm^^sr^^ respectively. IVermeii et al.l (l2002al) 
also measured these lines for several regions in 30 Dor with the 
LWS which has a beam 80". They found that the [O I] intensity 
is about twice the [C II] intensity. 

To estimate the contribution of [OI] in our 15' beam, we as- 
sume two co mponents within the b eam: 7 ) the smaller region as 
measured bv lPoglitsch et al.l (Il995h . and 2 ) more extended emis- 
sion. F or 1, we use the values as measured by 'Poglitsch et alj 
([19951) . For 2, we us e the differences between the line strengths 
as liste d by IVermeii et al I (|2002a) and listed by iPoghtsch et"aLl 
( Il995h . We assume that the ratio of [O I]/[C II] in the extended 
region holds throughout the 15' beam and scale that to the full 
[CII] measured in the large beam. This is a very conserva- 
tive estimate of the [O I] contribution because it probably over- 
estimates the [OI] line strength as the true ratio most likely 
decreases with distance. We thus estimate the contribution of 
[OI] 63^m to be ~20% of the total gas cooling rate. This num- 
ber should be lower in other regions in the LMC since 30 Dor 
represents the most extreme SF region. We conclude that, for 
our 15' beam, the [C II] is representative of the total gas cooling 
rate, but may slightly underestimate the cooling rate in the most 
extreme cases. 

Of the SF regions, 30 Dor has the lowest efficiency as it is 
the most intense SF region in the LMC (See Tab. |2]i. In con- 
trast the second brightest SF region. Nil, has a higher than aver- 
age efficiency. The regions defining 30 Dor and Nil were cho- 
sen by constructing a rectangle centred on the brightest pixel in 
the 24/im image. The edges of the rectangles approximate three 
times the mean S24 level of the diffuse regions. We have included 
the resulting rectangles in Fig.[T^. One puzzling result of this is 
that on average the diffuse regions exhibit a lower efficiency than 
do the SF region. It should be noted that most of the diffuse re- 
gions is faint and therefore affected more by the high noise in 
the BICE map. 

7.4. The limited sensitivity of tlie BICE map: validity of 
observed efficiency variations 

Before we can draw conclusions on the variations observed in 
the observed photo-electric efficiency parameter, we explored 
the possible existence of a systematic variation of the efficiency 
simply as a function of Go by running a set of simulations. The 
BICE map has limited sensitivity and as a result a large part of 
the diffuse medium remains undetected or has flux levels com- 
parable to the noise level (See Fig.|3]l. Therefore, extra care must 
be taken in deriving representative values of the [C II] intensity 
in the diffuse medium. In these simulations we assume that Stir 
is a proxy for Go. We simulate observed [C II] maps by applying 
a function of the form e=f(STiR)- The Stir map is constructed 
using the Spitzer images and the uncertainty in this map is neg- 
ligible compared to the BICE map. After this we add noise ac- 
cording to the noise level of the BICE map. Fig. Q compares 




Figure 7. Huxes as observed in the BICE maps versus those ob- 
tained from a simple simulation (section l7T4] i. The figure shows 
the observed histogram of [C II] surface-brightness (solid), the 
bin width is 1 cr in [CII]. The best fit model is shown in black 
(text for details). The dark grey line indicates the result of a sim- 
ulation with a moderate decrease in efficiency as a function of 
TIR. The light grey line corresponds to a model with an decrease 
by a factor 4, the latter simulation is clearly inconsistent with the 
observations as it over predicts the number the number of bright 
pixels. 



the observed histogram of [C II] values with several predicted 
histograms obtained following the simulation outlined above. 
We have simulated maps using functions of two forms linear) 
f(STiR) = a + 6*Stir and power-law) f(STiR) = a*S?^jf{. We de- 
termine the parameters which best approximate the distribution 
of the observed S[cii] distribution. Some results of these simu- 
lations are shown in Fig.|7] The closest match is obtained when 
assuming a constant efficiency across the full range of Go. The 
data do not exclude modestly higher efficiencies in the faintest 
regions by about a factor of two,. Steeper gradients are excluded, 
since they clearly over-predict the number of pixels with values 
between 3 and 5 <t in the BICE map. At a given value of Go 
we do find a significant spread of S[c ii]/Stir in our simulation. 
The spread, which is independent of Stir, has a magnitude of 
roughly a factor of two. The main conclusions of the simulations 
are as follows: 



- The value of S[cii]/Stir is independent of Stir over the 
range from 3 10~^ to 3 lQ~^erg s~^cm~^sr~^ . 

- The mean value of e across the LMC is 0.45 % 

- There is a spread around this value between 0.3 and 0.6 

- There is a modestly higher value of S[cii]/Stir, up to ^ 1% 
in the faintest regions 

- The only region exhibiting a systematically lower value is 
the 30 Dor region with a mean value of Spuj/STiR, of 0.35 
%. 
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Figure 9. Diffuse and SF data plotted over Fig. 17 from IStacey et alJ (|1991). .The solid line is a power law fit to our LMC data. 
Tracks of constant PE efficiency are also plotted (dashed lines). The first track on the left has an efficiency of 5% and the last track 
on the right has an efficiency of 5 • 10^"^%. The interval between each line is a step of 1.5. Colours and symbols as in Fig.H] see 
online version for colour 
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S^i^ [10"^ erg ^^cnr^sr] 

Figure 8. S[c nj plotted as a function of Stir for the diffuse and 
SF pixels. The solid line is a power law to the data. The dashed 
line is a linear fit to the data with the y intercept set to zero. 
Colours and symbols as in Fig. HI Colours and symbols as in 
Fig. El see online version for colour. 



8. [C II] Emission, Phiotoeiectric Efficiency and 
Radiation Fieid in thie LMC 

Here, we examine the variations in the observed [C II] emission 
and PE efficiency as functions of radiation field. 



8.1. S 



[CII] 



and Radiation Field 



We can take Stir as a proxy for the UV radiation field in the 
case that most of the power absorbed by the dust is in the UV, 
and the dust radiates isotropically. We plot Spii] as a function 
of Stir in Fig- 13 The [C II] emission increases as a function of 
Stir, and flattens at the highest radiation fields, which is clearly 
seen in the power-law fit. This flattening is dominated by the low 
efficiency observed in 30 Dor. We also tried to fit a straight line 
to the data with the y intercept set at zero. The data is clearly 
better described by the power law. The difference between the 
power-law and the linear fit at the highest values of Stir is 30%, 
which is more than we can comfortably explain by missing line 
emission from other cooling fine structure lines (Sec. l7.3l l. The 
power-law fit yields the following prescription for the [C II] sur- 
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face brightness as a function of total infrared surface brightness 
throughout the LMC: 

S[c II] = 1.25 -10-3 

where surface brightness is given in erg s~^cm^^sr^^. 
The equation is valid for Stir between 3.2 and 33 • 
10~''erg s~^cm~^sr^^, the range from which the fit was made. 

We interpret the flattening at the highest TIR in Fig. [8] as 
a decrease in the PE heating rate. In high radiation field and 
high temperature environments, grain charging effects become 
important, the PE heating efficiency decreases and thus the line 
cooling drops. This flattening is also ob served in KAO dat a 
of Galactic and extragalactic regions by IStacev et al.l (Il99ll) . 
In Fig. 17 of their paper, they plot Sjcuj as a function of the 
FIR surface brightness (which they call xfir) normalised to 
2 • s~^cm~^sr~^. We show our data, their data and 

Eq.|4]in Fig.|9l We have converted our TIR values into yfir by 
assum ing a factor of two for TIR to FIR as given bv lHunter et al.l 
(1200 ll) . We also plot tracks of constant PE efficiency in Fig. |9] 
(dashed fines) from e = 5% (left) to 5 • 10^'^% (right); where 
e = 5% is c hosen because it is close to the highest efficiency in 
flie theory of iBakes & TielensI (11994 . 

Our data thro ughout the LMC follow the trend of 
IStacev et al.l ( Il99ll) and extends to the lower left portion of the 
plot. To understand why our SF points do not occupy the up- 
per right portion of the graph, one must consider that our beam 
size (~ 200pc) undoubtedly entails considerable mixing of the 
phases of the ISM. Stacev et al. (1991) include measurements of 
30 Dor (indicated on Fig. |9]l. We use this to gauge how beam 
size affects this figure. The measurement of 30 Dor from our 
data is indicated with a star. 30 Dor in the BICE beam has sig- 
nificantly lower IR and [C II] surface brightnesses, by factors of 
about 125 and 15 respectively. Note, that because of the much 
smaller beam of t he KAO, and because the Galactic regions are 
closer, the data of IStacev et alj ( 1199 lb probe much smaller spa- 
tial scales. On these scales, the contribution to cooling from the 
[O I] line as compared to the [C I I] line might be im portant. If 
the [O I] line were included in the lStacev et al.l (11991) data, the 
points at the upper-right side of the diagram would move up. 
Even so, inclusion of these other lines would most likely not 
be sufficient to move the data points up to the efficiency of the 
LMC, as this would require them to be moved up the y-axis by 
at least an order of magnitude. Such cooling contributions from 
other lines, even in the most intense regions, are not expected. 
As with the decreasing trend of efficiency in the Stacey data, we 
note that for our data, 30 Dor shows the most prominent decrease 
in efficiency. Perhaps this reflects the fact that 30 Dor is so bright 
that it dominates the emission even in the large beam. 

One must note though, that the other variables which control 
the PE heating (i.e. gas temperature and electron density) are 
not constant throughout Figs. [8] and |9] Perhaps the considerable 
spread in [C II] at any given Stir in Fig.[8]and can be attributed 
to those other variables. 



9. Correlation Between Radiation Fieid and 
Electron Density 

In this section we use the data and the physics of PE heating to 
derive values for Gq and nc. We find that the observed constant 
efficiency translates into a tight correlation between the density 
of the radiation field and the density of electrons. We discuss the 
physical interpretation of this correlation and in particular we 




Figure 10. Go as a function of Uo for the SF pixels as estimated 
from Eqns. [3] and |5] The solid and dashed lines are fits of the 
equation Go = £,n^ ■ The solid line is fit with both ^ and p as free 
parameters and results in ^ = 2, 200 and p = 1.1. The dashed 
line is a fit with p set to 4/3 and ^ as a free variable. The result is 
i = 5,300. 
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Figure 11. Values of Go for each pixel derived from the dust 
temperature assuming graphitic grains for the dust composition 
versus values of Go given by Eq. |5] The solid line denotes 
agreement between the two methods. Colours and symbols as 
in Fig. m see online version for colour. 

compare our results w ith those of I Young Owl et al] (|2002|) and 
lMalhotraetan(l2001h . 

To first order the observed efficiency is constant throughout 
the LMC. This constancy of 0.45% translates into a constant 
value of 7 sa 4 10* (Eq. A typical value for PDR gas tem- 
perature is T '--^ 300K. To acc ount for the co ntribution from the 
diffuse ISM (T - 50 - 100K lTielensll2005 ')) we adopt a value 
of T = 75 K. Assuming an average temperature of 75 K across 
the SF pixels we find that all of the LMC regions, averaged over 
our 225 pc regions here, are typified by Go/ne « 5000. 

Several explanations for the constancy of the observed effi- 
ciency and, as a consequence, the derived interrelation between 
Go and ne come to mind. 

i) If, in our large beam, the radiation field is dominated by 
similar PDR regions, then the number of PDRs in that beam 
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Gg [habing] 

Figure 12. Comparison of the interrelation between Go and no 
that we find and that is found within galactic PDR regions by 
Youn g Owl et al. (2002). The diamonds indicate the data from 
Young Owl et al.l (12002 ). while the dashed line is their best fit. 
The drawn line shows ng as a function of Gq for each fixel in 
our analysis. Here no is derived from ne (Fig.fTOb assuming that 
the electron density equals the carbon density. 



will determine variations in the total radiation field, but will not 
change the intrinsic PE efficiency. 

ii) If, in our large beam, the light is dominated by emis- 
sion from the diffuse medium, PE efficiency should remain fairly 
constant as the physical conditions within this medium are rela- 
tively invariant. 

Hi) The value of Gol^e is constant throughout the differ- 
ent types of media tha t make up the LMC. In fact, for several 
Galactic PDR regions. Young Owl et alj (12002! ) do indeed find 
that Go and ng scale with each other They studied the PDRs of 
a sample of reflection nebulae. Their observations of FIR atomic 
fine structure lines and the FIR continuum allowed them to ob- 
tain estimates of radiation fields and gas densities (no). 

It is interestin g to compare our r esults in m ore detail with 
those obtained bv lYoung Owl et alj (l2002h and iMalhotra et alj 
(I2001h . The latter measured line emission ratios across entire 
galaxies for a large sample of galaxies with varying morpholo- 
gies. They used PDR models to obtain values for radiation fields 
and gas densities. Both then fitted power laws to their data . For 
a function of the form. Go = f 'ng.'Young Owl et al l (|2002|) find 



= 0.09 and p = 1.2, and IMalhotra et alj (1200 Ih £' = 0.2 3 
and p = 1.33 (estimated from Fig. 1 1 of Malhotra et al.1 ( 1200 lb ). 
In other words, they also find that the density scales roughly lin- 
early with the radiation density. 

Here, we estimate Gq and for independent pixels in the 
LC. We do this only for the SF-regions. One, because these 
are the brightest, reliably detected regions and two, because the 
conversion from Stir to Go uses the assumption of central il- 
lumination which is more likely to hold in those regions. We 
use the observed efficiency to invert the efficiency equation of 
[Bakes & Tielens ( 1994, Eq. (|3]l). This yields values of 7 for each 
pixel. We measure the Go using Stir, and assume a reasonable 
value for T. Thus we can derive values of ng for each pixel. 



Following t he exa mpl es of iMeixner et al.l (|1992|) ; 
l^ung Owl etiD (I2002h : ISteiman-Cameron et all (Il997h 
we use the measured TIR surfaces brightnesses and assume a 
certain geometry in order to estimate values of Gq. Assuming 1 ) 
that the illuminating sources are at the centre of each pixel, 2) 
that the dust resides at edges of the pixels, 3) 100% conversion 
of FUV radiation to IR radiation, 4) and that the dust radiates 
isotropically, then 



Go — G; 



s^G S'tir 



(5) 



with Cs-,G = 4 10"**, Go in units of the one dimensional habing 
field (1.610^'^ ergss^^cm^^) and 5tir in ergss^^cm^^. 
Other geometries; like the dust residing halfway between the 
source and the edge, homogeneous distribution of sources and 
dust or uniformly illuminated dust; changes Cs^g only by a 
small amount. The Go thus calculated is shown in Fig. [10] ver- 
sus the derived electron density. We find a best fit power-law of 
the form Go\ habinq] = 2 20 0ne^ -^[cm~^]. Note, that the re- 
sults of Malhotra et al.1 ( 1200 lb and lYoung Owl et all (l2002h are 
expressed in no while we use ng and thus they differ by the ioni- 
sation fraction. 

In order to numerically compare our result we convert ng 
to no assuming that the average electron density is dominated 
by diffuse-ish, outside of the H II regions, regions. In this case 
the prime donor of free electrons is carbon which is the most 
abundant element with an ionisation potential (11.3 eV) below 
that of hydrogen (13.6 eV). Taking into account the metallicity 
of the LMC we find a conversion factor (H/C) from ng to no of 
6000. It can b e seen in Fig. [T2l t hat ou r values follow the trend 
est abhshed bylYoung Owl et al.1 ( l2002h well. 

lYoung Owl et al.l (|2002|) and iMalhotra et all ([2001') pro- 
pose that the correlation can be physically interpreted as a 
bala nce of gas pressures be tween the PDR and the HII re- 
gion. l^un^OwTeFaridlOO^) present a simple analytical model 
assuming this balance of pressures and they find that Go 

should scale w i th n^ ^. For both lYoung Owl et al.l (|2002|) and 
iMalhotra et all (I2OOII) . the power of 4/3 fits to within the un- 
certainties of their data and fits within the envelope of the un- 
certainty of our data (Fig. [TOl i. We have performed an F-test to 
determine whether the fit with p as a free parameter is an im- 
provement compared to strict pressure balance (See Fig. [TO] i. 
The shallower power-law (p=l.l) does decrease the reduced 
significantly although it remains well above unity. 

The cause behind the correlation is unclear, especially one 
that would hold over such a range of scales. It is difficult to en- 
visage pressure-balance between the powering H II regions and 
the 200 pc regions of space that we are sampling. This is even 
more true for the results of Malhotra et al. (2001), who sample 
entire galaxies and still find a strong correlation between the av- 
erage radiation field and the average density. Such a pressure 
balance would work if each single region or galaxy is dominated 
by a single (or at most e few) cluster(s) of young stars that cause 
a single prominent PDR to plough into the containing molecular 
cloud. 



9. 1. Dust Temperature in the LMC 

In the previous sections, we used the measured TIR surface 
brightness as a proxy for radiation field. There could be a worry, 
though, that Stir can vary without actual variations in the radi- 
ation field. This may be due to differences in the amount of and 
characteristics — emitting material along the line of sight as a re- 
sult of varying densities or a varying thickness of the LMC. We 
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Table 3. Results of Efficiency Calculations 



Region 


es 


£24 


£70 


£160 


x' 


Integrated over LMC 


0.02 


0.00 


0.02 


0.01 


1.1 


SF regions 


0.06 


0.00 


0.00 


0.01 


2.0 


30 Dor 


0.05 


0.00 


0.00 


0.01 


4.7 


Nil 


0.05 


0.04 


0.00 


0.00 


2.8 


Diffuse regions 


0.00 


0.11 


0.02 


0.02 


0.8 



explore these concerns by comparing Go as derived from Stir 
to Go as derived by a rough calculation of the dust temperature 
as measured by the 70 /im to 160 /im ratio. This ratio in not 
affected by varying amounts of material along the line-of-sight. 

We derive a dust temperature indicator assuming modified 
black-body radiation with an emissivity of the following func- 
tional form: Q oc , with j3 = 1.5. This is a crude approxima- 
tion, since in reality dust grains do not emit as modified black- 
bodies. The dust temperature will vary with grain-size and also 
with grain composition. Moreover, a part of the emission we ob- 
serve at 70 /im is not due to dust at an equilibrium temperature 
(See Sec. |7]), but arises from stochastically heated grains. For 
these reasons, the absolute value of the derived temperature is 
not very reliable, but the relative temperatures are. In spite of 
the above caveats, the temperatures that we find are all-together 
reasonable and compare well with the results of more in-depth 
studies. The mean dust temperature we find (Td = 24 .0 ± 2.2K) 
is consistent with the ones found bv lSakon et al.l (12006.) (23.9 isT) 
andlAguirreetaD (l2003|) (25.0 ± 1.8 K). 

Us ing the relationship between Go and given by iTielensI 
(l2005b for graphite grains, we calculate Go(Td) for each pixel. 
We plot Go(Stir) versus Go(Td) for each pixel in Fig.[TT] The 
solid line denotes perfect correspondence. The correspondence 
between the two Go indicators is clear and thus we conclude that 
Stir is indeed a valid proxy for Gq. 



10. Photoelectric Efficiency and Grain 
Distribution in thie LMC 

Theoretical studies have shown that the extent of grain charg- 
ing and therefore the efficiency of PE heating is not only de- 
pendent on environmental conditions, but also on grain size and 
grain species. Interstellar molecules such as PAHs and the small- 
est gr ains contribute most extensively to the PE heating (Watson 
119721; [Bakes & Tielens 1994; W eingartner & Drainell2001 ). For 
example, .Bakes & Tielens(.1994i) have found that approximately 
half of the heating is from grains with less than about 1 500 C- 
atoms (^ 15 A). 

Qualitatively, the scaling of efficiency with grain size is due 
to the fact that the ionisation rate is approximately oc Nc since 
the FUV absorption cross section is approximately cx Nc. But, 
the recombination rate increases more slowly with the number 

1/3 

of carbon atoms as it scales with the grain size (about NJ ) 
. Therefore, the fraction of ionised grains increases with grain 

2/3 

size, scaling by about N J . 

There has been a plethora of theoretical studies on grain 
species and PE heating, but relatively few observationally based 
studies. Therefore, we now examine PE efficiency as a function 
of grain species by using each band as a tracer of grain abun- 
dance. 



10.0.1 . A Calculation of Efficiencies 

Authors usually take [C II]/FIR as a general efficiency, encom- 
passing the contribution to PE heating of all constituent grain 
populations. But, as each grain component is expected to have 
different intrinsic PE efficiencies, it is interesting to isolate the 
contribution of PE heating from the various species. Therefore, 
to quantify the importance of the various g rain species to the P E 
heating process, we follow the example of Hab art et al.l (1200 ll) . 

Using IRAS data, Habart et al. (2001) quantified the amount 
of emission attributed to the PAH, VSG and BG populations in 
the p Opiuchi complex. Using ISO observations of the [CII], 
[Ol] and H2 lines, they then calculated epAH, evsG^ cbg with 
a linear combination fit of their gas cooling rate to their grain 
emission rates (further detail given below). They found that the 
PAH population is attributed with the highest photoelectric ef- 
ficiency, while the BG population is attributed with the lowest. 
They found (epAH = 3%, eysG = 1%, esc = 0.1%). 

To perform this calculation, we start with the definition of 
efficiency for a certain grain population,], which is given as 

^ rpE(pop.i) ^ S[cn](popj) 

Pabs(popj) Soinittod(popj) 

The substitutions of FpE and Pabs with S[c n] and Somittod were 
already discussed in Sec. 17.31 To solve for the total amount of 
emitted [CII], we multiply both sides of Eq.|6]with the denom- 
inator on the right side of the equation, and then sum over all 
grain populations: 

^[CII] = 5Z £pop jScmittcd(pOp j). (7) 
all pop 

Finally, if we consider the emission in the 8, 24, 70 and 160/im 
bands as arising from distinct grain populations, and that these 
are the main populations that contribute to PE heating, we ob- 
tain: 

S[CII] — fsSs + £24824 + £70870 + £l6oSl60- (8) 

The variable, £§, should reflect the of PAH efficiency, £24 of the 
VSG efficiency, £70 of the a combination of the VSG and BG 
efficiency and eigo of just the BG efficiency. 

To solve for the efficiencies of Eq. [8] we fit a linear combi- 
nation of the IRAC and MIPS data to the [C II] map with a 
minimisation. Further, we force the efficiencies to be between 
and 1 . We perform this calculation across the whole galaxy, for 
the 30 Dor and Nil regions and for the diffuse and SF regions 
in the LMC. All calculated efficiencies, along with the value 
for each fit are given in Table [3] 

Table[3]shows that with the exception of Nl 1, £24 and £70 are 
always zero and £§ and £i6o are non-zero. We show the values 
derived for the diffuse medium as well. However, these values 
should be taken with extreme care, because of large contribution 
of noise to these pixels. We have experimented with deriving 
grain-efficiencies also using simulated [C II] maps and find that 
the 2% value for the BG is found persistently but that the £-VSG 
grain is very sensitive to the exact noise characteristics of the 
[C II] map and can not be trusted. 

The values in Table [3] show that £3 is greater than figo for 
every region considered. This is quantitative proof that the PAH 
emission is more spatially correlated with the [CII] emission 
than the BG emission. This supports the interpretation put forth 
in the previous section, that the PAH population dominates the 
PE efficiency and plays an important role in the Photoelectric 
heating of the gas. 



Rubin et al: [C II] and dust emission in the LMC 



15 



The results for Nil obviously differ from the other regions 
considered as eigo — and as and £24 7^ 0. This might reflect 
offsets between the [C II] emission and the Spitzer bands due 
to the distinctive asymmetry towards Nil already discussed in 
Sec.E] 

11. Conclusion 

Using the IVIIR to FIR SAGE maps, and the BICE [C II] map of 
the LMC, we have, for the first time, conducted an observational 
study of PE heating and [C II] cooling in relation to spatially 
resolved grain emission throughout the LMC. 

Integrated throughout the entire LMC, the [C II] line ac- 
counts for 0.64 ± 0.01% of the total infrared (- 1.2% of the 
FIR). Applying a correction for the pixels below 2cr[c nj , we find 
that the [CII] fine accounts for 1.32 ± 0.01% of the FIR. This 
value is greater than that of normal and gas rich galaxies (with 
values typically from 0.1-1%), as found in other low metallicity 
galaxies. 

Distinguishing environments by Ha surface brightnesses and 
by location, we find that the [C II] line contributes significantly 
less to the TIR emission in SF regions versus diffuse ISM re- 
gions: 

- L[c iij/Ltir = 0.57 ± 0.01% for the SF regions 

- L[cii]/Ltir = 0.42 ± 0.01% for 30 Doradus 

- L[c II] /Ltir = 0.56 ± 0.03% for Nl 1 

- L[cji]/Ltir = 0.79 ± 0.01% for the diffuse regions. 

We also calculate the contribution of the total output of [C II] 
emission from the LMC from the same regions. We find that, al- 
though the SF regions have the highest surface brightness values, 
most of the LMC's [C II] emission originates from the diffuse 
medium: 

- L[cii]/L[cii],LMC = 31.5 ± 0.1% for the SF regions 

- L[cii]/L[cii],LMC = 5.58 ± 0.1% for 30 Doradus 

- L[cii]/L[cii]!lmc = 1.53 ±0.1%for Nil 

- L[cii]/L[cii],LMC = 68.5 ± 0.1% for the diffuse regions. 

To examine variations in PE efficiency within the LMC, we 
use S[cii] as a proxy for the total PE heating rate. We estimate 
that this assumption is valid for our 15' beam for all but the most 
active SF regions which might underestimate the total heating 
rate by at most 20%. We study how PE efficiency varies with 
environment using the observed values of Stir as an indicator 
of the local radiation field. As a function of Stir, S[cu] flattens 
at the highest value of Stir. The flattening trend is interpreted as 
a decrease in the PE efficiency for the most illuminated regions 
in the LMC. We provide a prescription for S[c n] as a function of 
Stir (Eq.|4|i by fitting a power law to the data. Such a decrease in 
efficiency is theoretically expected due to grain charging effects. 

Previous studies have found a correlation between these two 
parameters, with Go oc no^"^. This relation has been explained 
by invoking a simple model assuming pressure-balance between 
H II regions and the adjacent PDRs. Theoretically, the PE effi- 
ciency depends strongly on the recombination-rate, and thus on 
the ratio of Go/no. We thus calculate values for Go and us- 
ing the observed efficiencies. We convert the observed Stir to 
Go, assuming illumination by a central source in each 45 x45 pc 
pixel. We find that a similar scaling-relation between Go and iic 
holds for the LMC. It is unclear why the Stromgren sphere argu- 
ment should hold on the large ^^45 pc scale that we probe. 

We analyse observed PE efficiencies in relation to the grain 
component emission from each Spitzer band. We note that this is 



the first such analysis utilising spatially resolved grain emission 
components throughout an entire galaxy. From the correlations 
between observed efficiency and the grain component emission, 
and from a calculation of the PE efficiencies for each population, 
we show that the PAH emission is the most spatially correlated 
with the PE heating rate. We therefore conclude that it is the 
PAH population that dominates the PE heating process. 

The efficiency of the PE heating process is dependent on 
both environmental conditions, such as radiation field, and grain 
abundances. Our study has examined PE efficiency within the 
LMC without fully disentangling the extent of PE heating due to 
existence of grain populations favourable to PE heating and the 
extent due to environmental conditions favourable to PE heating. 
Disentangling the effects of both on the observed PE efficiency, 
however, is difficult as the regions where PAH populations are 
destroyed naturally have intense radiation fields which also sup- 
press the extent of PE heating. To break this degeneracy we will 
undertake detailed SED modelling of the different regions in the 
LMC in an upcoming paper. Using the SED models we can inde- 
pendently solve for radiation field values and PAH abundances. 
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